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Abstract 

In present paper we propose an approach based on examination of 
the structure of the general solution of equations of the type dy/dx — 
P{x,y)/Q{x,y), with P and Q polynomials only in y. Under the term 
structure we mean the dependency character of solution from arbitrary 
constant. We describe a common form of the structures for foregoing 
equations. In such a way one can obtain a differential-algebraic polynomial 
system for undetermined parameters of the structures. The successful 
solution of this system automatically leads to finding the general solution 
of ODE's. We demonstrate on examples that proposed method gives, as 
a first step, the systematic way for obtaining new integrability conditions 
and general solutions for various families of rational first-order ODE's. 



1 Introduction 

Equations of the type dy/dx = P{x,y)/Q{x,y), with P and Q polynomials in 
y arise in many different practical fields. There are some methods for finding 
the solutions for such type of equations. First of all it is the famous method of 
integrating factor. A remarkable method based on the knowledge of the general 
structure of the integrating factor was developed by Prelle and Singer for the 
case when P and Q are polynomials in both arguments x and y (see also Q ) . 

In present paper we propose an approach based on examination of the struc- 
ture of the general solution of equations with P and Q polynomials only in y. 
Under the term structure we mean the dependency character of solution from 
arbitrary constant. In theory of differential equations (especially linear) the use 
of solution structure is very fruitful. Such type of approaches allows us, for 
example, to obtain the general solution from the knowledge of definite number 
of particular solutions (superposition) or to use the method of undetermined 
parameters and so on. 

Below we will see that many of practically important equations have suf- 
ficiently easy structures. The knowledge of the structure allows us to put up 
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different kinds of solving strategies. In present paper we expose the preliminary 
results, which are directed on demonstration of the method abilities. 

The paper is organized as follows: in section 2 we consider the base of the 
approach and discuss some possible solving strategies. In section 3 the well- 
known example of structure of Riccati equation is discussed from standpoint of 
proposed method. Sections 4-6 are devoted to specific examples, which expose 
some method features. 



2 Base of the approach 

The general solution of the first-order ODE 

(1) 

is the function y{x, c) of two variables x and c, where c is an "arbitrary constant" . 

Let is given a function of two variables y{x,c). Is there an ODE of type 
which general solution is given y{x,c)7 To answer this question it is convenient 
to use well-known connection between first-order ODE's and linear first-order 
PDE's. If C{x,y) is a solution of linear PDE 

^^ + /(x,y)^^=0, (2) 

then the general solution (in implicit form) of the equation (|l|) is 

C{x,y)=c, (3) 

and it being known that 

y{x, C) = c. (4) 

So, if for given function y{x, c) we solve algebraical equation in respect to C 
then from (|^) (under obvious restrictions) we obtain that the equation (|l|) with 

9y 

has as the general solution exactly the given function y{x,c). 
The examination of the expression (|5|) shows that if 

C{x,y)^F[— - + aln— -+/?arctan— -, (6) 

qi{x,y) q2{x,y) qz{x,y) 

where F is an arbitrary differentiable function; pi, q.i (i — 1, ... ,3) are poly- 
nomials in y; a and (3 are constants, then f{x,y) is rational function in y. If 
we seek solutions for equations with complex f{x,y), then presence of arctan 
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term is needless. The corresponding expression for f(x,y) here is as follows 
(a = /? = 1 for short): 

f{x,y) = - [qi{x,yf-^^^^^^q2{x,y)qi{x,yf + 

2dp2{x,y) . ^ . x2 
+ qi{x,y) q2ix,y)p3{x,y) - 

- gi {x, yfp2{x, y) ^'^^^2^ '^^ qaix, yf - 

- qi{x, yfp2{x, y)S^^p,{x, yf + 

+ ^^^^r^q2{x,y)p2{x,y)qi{x,y)q3{x,yf + 
+ ^2 g2 (a:, y)p2 {x, y)qi {x, y)p3 {x, yf - 

-Pi{x, y) g2 {x, y)p2 {x, y)q3 ix,yf - 

ax 

- Pi{x, y)~^^^~1'^i^' y)P2{x, y)p3{x, yf + 
+ q2{x, y)p2{x, y)qi (x, yf^^^^^q3{x, y) - 

- q2{x, y)p2{x, y)qi{x, yfpsix, y)^^^^]/ 

2dp2{x,y) 2 , 

m{x,y) — — q2{x,y)q3{x,y) + 

2dp2{x,y) / X / ^2 

+ qi(x,yj 1 q2(x,y}p3(x,y) - 

dy 

2 / ^dq2{x,y) 2 

- qi{x,y) P2{x,y) — — — q3{x,y) - 

- qiix,yfp2ix,y)^^^^^j^P3{x,yf + 

dy 

+ ^^^i^' q2{x, y)p2{x, y)qi{x, y)q3{x, yf + 
dy 

+ ^^^^-^r^q2ix,y)p2{x,y)qi(x,y)p3{x,yf - 
dy 

- P\{x, y)^^^^^^^q2{x, y)p2{x, y)q3(x, yf - 

dy 

- Piix, y) '^^ q2{x, y)p2(x, y)p3{x, yf + 

+ q2{x, y)P2{x, y)qi{x, yf-^^^^^q3{x, y) - 

dy 

I \ I \ I \2 I \dq3(x,y) . 

- q2[x,y)p2{x,y)qi{x,y) P3{x,y) (7) 

dy 



We can interpret the relation from different points of view. If we, for 
example, specify definite expressions for pi and g^, we can obtain a family (de- 
pending on parameters of the polynomials) of first-order ODE's with known 
(predetermined) structure of general solution which in implicit form is given by 
. Such interpretation as a first step of analysis is very essential in methodical 
sense. 

Practically more important is an attempt to solve inverse problem: guided by 
some consideration of orders of polynomials we can try to equate given rational 
function f{x,y) with a right-hand of (0) with undetermined functions - coeffi- 
cients of polynomials. In such process when we equate polynomial coefficients 
we obtain a differential-algebraic polynomial system of ODE's for undetermined 
parameters. The successful solution of this system (we have to stress here that 
it is sufficient to find particular solutions) automatically leads to finding the 
general solution of ODE of type (P with given rational function f{x,y) as 
right-hand term. Since in intermediate stage a polynomial system arises, there 
is hope that we can handle (simplify) this system by existent CAS elimination 
algorithms, e.g., by using casesplit from ODEtools in Maple. 

There are other possible variants of solving strategies, as e.g., on that way it 
is easy to reveal the structure of integrating factor and further to seek solution 
by means of undetermined parameters. 

At last under known structure of solution (in such a way it is sufficient 
to ascertain the fact of existence of a solution of the system for undetermined 
parameters with given f{x,y)) we can try to find the necessary number of partic- 
ular solutions of initial ODE and further to evaluate undetermined parameters 
via the set of particular solutions (non-linear superposition). 

Let us count roughly the orders in y of numerator and denominator of ra- 
tional function f{x,y) in (|^). Let us denote orders of polynomials pi and qt as 
Up. and rig.; N = X]?=i('*pi +"-9!)- 1^ is obvious, that order of polynomial P in 
numerator 

np < max{(iV + Uq^ + n^., - rip^ - Up^), 

{N + Ug^ - flg^ - rip^ + Up^), 
{N + Ug^ - Up^), 

{N ~ng., +7ip3), 

(N + Ug, - Up,)} 

and order of polynomial Q in denominator 

nq <np - 1. 

The maximal possible number of parameters (the general number of coefficients 
of polynomials pi and qi) in the right-hand of (0) is + 6 (A^ -I- 4 when arctan 
term is absent, or A -I- 2 when In term is absent too). As every rational function 
may be reduced by a nonzero factor and by virtue of elementary identities of type 
\n{a{x)4>{x, y)) = In a{x) + In (l){x, y) then the real number of free parameters is 
only A^ -f 1. 
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As we would like to obtain an arbitrary rational function assigned in advance 
which nominally would have np + nq + 2 (or np + nq + 1 free) parameters 
we conclude that generally speaking there are a deficit of free parameters. It 
means first of all that for the most cases the classical classification of ODE's on 
orders of polynomials |^ does not coincide with classification of structures of 
solutions. To one classical class of ODE's (e.g., Abel equations) corresponds the 
set of structures, depending of parameters of specific equation. The mechanism 
of structure selection (determination of degrees for polynomials pi and qt) is 
similar to proposed in [^]. 

From another side, as a consequence of a deficit of free parameters, above 
mentioned system can be solved if some definite conditions (integrability condi- 
tions or constraint) for functions-parameters of given equation are held. 

In the next sections we demonstrate on examples that proposed method 
gives, as a first step, the systematic way for obtaining new integrability condi- 
tions and general solutions for various families of rational first-order ODE's. 



3 Riccati equations 

The structure of solution is known for only a few non- linear equations. The 
most prominent example is the Riccati equation, where 

, , bO(x) c — aO(x) , . 

2^^'^^~mH tH 8 

bl[x)c — al[x) 



or 

_ al{x)y + aO{x) 
'^'^ blix)y + bOix)- ^ ' 

If we substitute (^) into (^, we obtain that this structure corresponds to 
ODE's family with 

^ ^bl{x)~al{x)^ 
■'^"''^^ -al{x)bO{x) + bl{x)aO{x)^ 

i^fL5Mw(x) + ^^61(x)-al(x)^^-aO(2;)^^ ^ 



~al{x)bO{x) +bl{x)aO{x) 
^^bO{x) - aO{x)^^ 



(10) 



-al{x)bO{x) + bl{x)aO{x) 
The general form of the Riccati equations is 

^ = X2{x)y^ + Xl(x)y + XO(x) (11) 
dx 

If we now equate rig ht-hand sides of (^ and (|ll|) we obtain the following ODE's 
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system: 



. — r — 'H .r) — (-(1 .r — ; — - , ^ 

^-al(x)60(x)+61(x)aO(x) 



-al{x)hO{x) + hl{x)a{){x) 

^60(.T)-aO(x)^ 



-al(a;)feO(x) +61(x)aO(x) 



= Xl{x); 

= XQ{x).} (12) 



We have to solve this system for al(a;), aO(a;), and 60(x). Let us simplify 

the system with help of differential elimination procedure, e.g., by using casesplit 
from ODEtools in Maple (everywhere in this paper we use casesplit with the 
default rankings). As a result we obtain some cases, the most interesting of 
which here is 

[al(x) = (-aO(x)X2(x)60(x) - ^^^bl{x) + 

ax 

+ Xl{x)bl{x)aO{x) + aO{x)^^^)/{bl{x)XO{x)) (13) 

= (^-bl{x)XO{x)aO{x)^^^^bO{x) - 2bl{x)X0{x)^^^X2{x)b0{x) + 

dx-^ dx ax 

+ bl{x)XQ{x)aQ{x)'^^^ - bl{x)XQ{x)Xl{x)'^^^aQ{x) + 
dx/ dx 

+ bl{xfXQ{x)Xl{x)'^^^ + bl{xfXQ{x)^^^^aQ{x) - 

_ 2{^^^fx(){x)aO{x) + 2XQ{x)X2{x)bQ{x)aQ{x)^^^ + 
ax dx 

^dblix) daiXyx)^^, , ^dX{){x) ^dbl{x) 

+ W(a.)^^aO(x)X2(a.)60(a=) + 61(0.)^^^^ - 

ax dx dx 

- W(x)2^^^Xl(x)aO(x))/(XO(x)W(x)2), (14) 
dbO{x) 2 



= (-X0(x)61(x)^ + bl{x)bQ{x)Xl{x) 



-bQixfX2(x) + b(iix)^^^)/blix)\ (15) 

ax 

where [61 (x) ^ 0, XO(x) 7^ 0] (16) 

Here the number of undetermined parameters of the structure is superfluous, so 
there are not additional (integrability) conditions, except obvious restrictions 
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(p.6[). We can conclude that any Riccati equation has the structure in form 
of (||). To express the unknown parameters al{x), aO{x), bl{x) and 60{a;) via 
parameters of Riccati equation we have to solve among others almost the same 
equation (|l^) but here it is sufficient to find its particular solution. 

4 Abel A type equations 

Let now consider by the same way the following slightly more complicated struc- 
ture 

C(a;, y) = al{x) y + bl{x) + ln(a2(a;) y + b2{x)) (17) 

which corresponds to ODE's with 

, ,dal(x) , , ,da2(x) dbl(x) ^, , 
f{^, y) - - {—r^a2{x) + {—^ + —^a2{x) + 
dx dx dx 

dal(x) , db2{x) , dbljx) ^^^ 

{al{x)a2(x) y + al{x)b2(x) + a2{x)). (18) 
This is sub-family of Abel A class of equations (see 1^): 



dy_ _ X2{x) y'^ + Xl{x)y + XQ{x) 
d^ ~ Yl{x)y + YO{x) ■ 

Unknown parameters here are defined by the following system 



da2{x] 
dx 



dal{x) 
^ dx '^'("^ 


= --X2{x 


dbl(x) ^, , dal(x),^, , 
■+ rf. '^^^^^^ dx ''^^^ 


= ~Xl{x 


db2ix)^dblix)^^^^^ 
dx dx 


= -XO{x 


al{x)a2{x) 


^Ylix); 


al{x)b2{x) +a2{x) 


= YO{x)}. 



(19) 



(20) 

In contrast to the Riccati equations here solutions of the system not always 
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exist. One of the cases here is as follows 

[h2{x) = {-XO{xfYl{x)'^ + Yl{xfXO{x)Xl{x)YO{x) - 

- Yl{xfYO{x)XO{x)^:^^^ + Yl{xf^^^^YO{xfXl{x) + 

ax ax 

+ Yl{xfYO{xfXQ{x)^^^^ - Yl{x)^^^^YO{xfXl{x) - 

ax ax 

- Yl{x)Xl{x)X2{x)Y0{xf - Yl{x)^:^^^^X2{x)Y0{xf + 
+ X2{xfY0{x)^ + ^^^^X2(x)y0(x)4)/ 

{-Wl{xf^^^^X2{x)Y0{x) + Yl{x){^^^^fYO{xf + 
ax ax 

+ AYl{x)X2{xfY0{xf - AYl{xfX2{x)Y0{x)Xl{x) + 
+ YlixfXlixf - 2Ylixr'^'-XmY0ix) + 

+ 4Yl{x)X2{x)Y0{xf^^^ - 2Yl{xf^^^Y0{x)Xl{x) + 

+ Yl{xf{^^^f + 2Yl{xf'^^^Xl{x)); (21) 
ax ax 



ax ax ax 

+ Yl{x)Xl{xf - 2Yl{x)X{){x)X2{x) - Xl{x)Y{){x)X2{x) - 

-yO(x)X2(,x)^)/ 
ax 

{Yl{xfX{){x) + YQ{xfX2{x) - Yl{x)Xl{x)YQ{x)); (22) 



a2{x) = {-Yl{xfXQ{x) - YQ{xfX2{x) + Yl{x)Xl{x)YQ{x)) / 

+ Ylix)'^^^^ + Xl{x)Ylix) - 2X2(x)YQ(x)); (23) 
ax ax 



al{x) = {2X2{x)Yl{x)Y0{x) +Yl{x)^^^^YO{x) - 

- Y\{xfX\{x) - 

ax 

{Y\{xfX^{x) + Y'^'yxfX2{x) - Y\{x)X\{x)Y^{x)y, (24) 
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dPYQjx) 
dx"^ 

+ 3Yl{x)Y0{x)^X2{x)^^:^^ 

ax 



dx 



Yl{x)YQ{x 



y,dX2{x) dYl{x) 



dx 



dx 



+ ri(.)3^^inMro(.) 



Yl{xfXQ{x)X2{x) '^^^}'''^ - Yl{xfXQ{x) 
dx 



dx dx 
dYQ{x) dYl{x) 



dx 



dx 



dXl{x) 



dx 



y^^^^.dXlixl dYKx, 
da; 



dx 



- Yl{xfY0{x)^X2{x) 



+ Yl{xfXO{x){^^^fYO{x) - Yl{xfYO{x) 
dx 



dx 

^dX2{x) dYQ{x) 



dx 

^2.^n..^ ^n>^^t'^X^f _ 2X2{xfYQ{xf + 
dx 



dx 



2Yl{xfYQ{x)X2{x){- 



+ X2{x)X{){x)Xl{x)Yl{xf + ?,X2{xfY{){xfXl{x)Yl{x) - 

- Yl{xfX2{x)Xl{xfYQ{x) - 2Yl{xfX2{xfXQ{x)YQ{x) - 

, . . . ,dYQ{x) dYllx) , , , , ^dYllx) 

- Yl{xf Xl{x)YQ{x) -r-^ + Yl{xf XQ{x)X2(x) -^YQ(x 

dx 



dx 



2,Yl(xfXl{x)YQ{x)X2{x)'^^^^ + ■iYl{x)YQ{xf X2{x] 



dx 

dYO{x) dYl{x) 



2Yl{xfXid{x) 



dX2(x) 
dx 



dx 

YQ{x) + 2Yl{xf^-^ 
dx 



dx 



dx 



X2{x)YQ{x) 



Yl{x)Y0{xfX2{x)'^^;^ + iYl{x)YQ{xfX2{x)Xl{x)'^^^^''^ 



dx^ 



dx 



d^Yl{a 



Yl{xfXO{x) -i^YO{x) + Yl{x)^Xl{x)YO{x)^ 

dx'' 



d^Yl{x) 
dx'^ 



^^^o.^.^dYlix),,_^^^^^^,^^^^^,dYl{x) 



Y0{x)'^X2{x){- 



dx 



+ Yl{xfXl{xr^^^ + Yl{xfXl{x)C' "^-^^ ^2 



Yl{x 



dx 

^dYO{x) dXO{x) 



dx 
dYO{x) 



dx 



Xl{x)+Yl{xfXQ{x)'^^^^)/ 



dx dx dx dx 

{-Yl{xfXO{x) - Y\{xfY{){xfX2{x) + Yl{xf Xl{x)Y{){x))]] (25) 



where [Yl{x)^^^^^ + Xl{x)Yl{x) - 2YQ{x)X2{x) - YQ{x)^^^^^ ^ 0, 



dx 

Yl{x) ^ 0, 

Yl{xfXQ(x) + YQ{xfX2{x) - Yl{x)Xl{x)YQ{x) ^ 0]. 



dx 



(26) 



So we have to consider the equation (25) as an integrahility condition for this 
case. If for some Abel A type equation the integrabihty condition (2^) is held 
true, then the undetermined parameters of the solution structure (17) are easily 



calculated (in quadratures!) from (|21|)-(24) and the general solution immedi- 
ately follows from (H). Note that in the case of Abel A type of equation we have 



five parameters, four of them here are free and fifth is defined by the integrabihty 
condition. 

The existent Maple ODEtools implementations do not find the solutions for 
equations considered here and below. 



5 Unclassified example 1 

With the following structure 

^ a2{x)y^ + al{x)y + aQ{x) 

b2{x)y^+bl{x)y + b0{x) ^ ' 

the more wealthy family of ODE's is connected: 

/(., y) = - a2ix/^) + (-a2(.)^ - + 

+ ^61(0.) + ^62(0=)) y' + {^bO{x) + '^b2{x) - 
dx dx dx dx 

.dbl{x) dal{x),^, , „, .dbQix) „, .db2(x). , 

dx dx dx dx 

, ,dbO(x) , dbl(x) dal(x), , , daO(x), , 

((al(x)62(x) - a2(x)61(x)) y^ + 

+ (2o0(x)62(x) - 2a2(x)60(x)) y + aO(x)61(x) - al(x)60(x)). (28) 

In general it is unclassified equations of the following type 

dy_ ^ XMx) + Xi{x) y^ + X2(x) + Xljx) y + XOjx) 

dx Y2{x)y'^ + Yl{x)y + Y0{x) ' ^ ' 
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Unknown parameters here are defined by the following system 

{^^62(.)-a2(.)^=X4(x); 
ax ax 

^ dbljx) db2{x) da2ixl .^dalix^ 

-a2{x) — ; alia;) — ; \ ; blix) H ; o2{x) — a3(x); 

ax ax ax ax 

da2(x),^, , daO(x),^, , .dbl(x) dal(x) , ^ , , 

ax dx dx dx 

,dbO(x) ^, ^db2(x) , 
-a2{x) — - aQ{x) — = X2{x)] 
dx dx 

_ dbOix) ^. .dbl(x) dal(x),^. . daO(x)^^, . , . 

~al(x) — - aO(x) — p-^ + — -^bO(x) H -^bl(x) = Xl(x); 

dx dx dx dx 

^, .dbO(x) daO(x),^, , s 
dx dx 

al{x)b2{x) - a2{x)bl{x) = Y2{x); 
2a0{x)b2{x) - 2a2{x)b0{x) = Ylix); 

aO{x)blix) - al{x)bO{x) ^ YO{x)}. (30) 



Here we examine one of simple cases which corresponds to the equation ( p9[ ) 
with X4:(x) — 0. Another more complicated case we will consider in the next 
section. 

[al{x) = Y2{x)/b2{x), (31) 
bl{x) = {2bO{x)Y2{x) + 2b2{x)Y0{x))/Yl{x), (32) 
aO{x) =Yl{x)/2b2{x), (33) 

dbO{x) 



dx 



= -l/2{4:X0{x)b0{x)Y2{x) + 4:X0{x)Y0{x)b2{x) - 

- bO{x)Yl{x)Xl{x) - bO{x)Yl{x)^^^^)/{Yl{x)YO{x)), (34) 

dx 



^^^^ = {-4Y2{x)b2{x)Y0{x)Xl{x) + AY2{x)b2{x)Yl{x)X0{x) + 
dx 

+ AY2ix)b2{x)Y0{x)^:^^^ - b2{x)Yl(x)^^:^^ - 
dx dx 

- 8X3{x)b2{x)Y0{x)^ - b2{x)Yl{x)^Xl{x) + ib2{x)Yl{x)Y0{x)X2{x))/ 
{-2Yl{x)'^Y0{x) + 8Y0{x)^Y2{x)) (35) 



a2{x) = 0, (36) 
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dY2(x) 

^ ' = {-Yl{x)^Y0{x)X3{x) -Yl{x)^Xl{x)Y2{x) 



dx 



- Y\{xY -^Y2{x) + \Y\{x)Y2(xYX'i^(x) + &\{x)YQ{x)Y2{x)X2(x) - 

dx 

- 4Y0{x)Y2{xfXl{x) + AY0{x)Y2{xf^:^^^ - AY2ix)Y0{x)^X3{x))/ 

dx 

{4Y0{xfY2{x) - Yl{xfYO{x)), (37) 
: {^8X0{x)Y0{x)Y2{xf + AY2{x)X0{x)Yl{x)^ + 



dx 



AY2{x) — --^Yl{x)YO{x) - Yl{xyXl{x) - Yl{x) 



dx dx 
+ 2Y0{x)Yl{xfX2{x) - 4X3{x)Yl{x)Y0{xf)/ 

{AY0{x)^Y2{x) ~ Yl{x)^YO{x)), (38) 
X4{x) = 0], (39) 

where [b2{x) ^ 0, YO{x) ^ 0, Yl{x) ^ 0, 4Y0{x)Y2{x) - Yl{x)^ ^ 0] (40) 

Here we have three equations (|37|)-(|39|) as integrabihty conditions. The un- 
determined parameters of the solution structure (Ujh are calculated (in quadra- 



tures too) from (31 )-(p6[) and the general solution follows from (^. Note that 
for this type of equations we have seven parameters, five of them here are free. 

Of course, this type of equation is defined by the structure with a2{x) — 
(or b2{x) = - see the first equation of (|30|)) but here we see that more simple 
cases are particular cases of more complicated structures. The main problem 
on this way is the limit of complication of structures, which is tackleable by 
existent differential elimination procedures. 



6 Unclassified example 2 

Examples in sections 4,5 lead to general solutions in quadratures. In some 
tackleable cases we can receive solutions with some special functions. As such 
example let us consider one of the cases of simplification of system for structure 
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( p7| ) of previous section. The case is as follows: 

[al{x) - (Yl{xfa2(x)^^^Y2{x) - 2Yl{xfY0{x)Y2{x)'^''^^'' 



dx dx 
+ Yl{x)^a2{x)Xl{x)Y2{x) - 4Yl{x)a2{x)X0{x)Y2{x)'^ - 

~ 'iY0{xya2{x)X4:{x)Yl{x) ~ 4:Yl{x)Y2{x)Y0{x)a2{x)X2{x) + 



+ 8Y2{x)^Y0{x)^^^^ + AY2{x)^Y0{x)a2{x)Xl{x) + 
dx 



+ 8Y2{x)a2{x)X3{x)Y0{x)^ ~ 4Y2{xfY0{x)a2{x /^^^''^ )/ 

dx 

{8Y0{xfX4:{x)Y2{x) - 2Y0{x)X4{x)Yl{xf), (41) 



bl{x) = i'iY2{x)'^Y0{x)a2{x)b2{x)Xl{x) +Yl{x)'^Y2{x)b2{x)a2{x)Xl{x) + 
+ Yl{xfY2{x)a2{x)b2{x)^^^^^^ + 2Y2{x)Y0{x)X4{x)Yl{xf - 



2Yl{x)^Y2{x)Y0{x)b2{x)^^^ + 8Y2{xfY0{x)^b2{x)^"'^^'''^ 



dx dx 

- 'iYl{x)Y0{x)^XA{x)a2{x)b2{x) - 'iYl{x)b2{x)Y2{x)^a2{x)X0{x) - 

-4Yl{x)Y2{x)Y0{x)b2{x)a2{x)X2{x) - 8Y2{xfY0{xfX4:{x) - 

dYO(x) 

- 4Y2{xfY0{x)a2{x)b2{x) + 8Y2{x)a2{x)b2{x)X3{x)Y0{xf)/ 

dx 

{-2Y0{x)a2{x)X4:{x)Yl{xf + 8Y0{xfa2{x)X4:{x)Y2{x)), (42) 



aO{x) = (Fl(x)3a2(a;)^^^ + 8Y2{x)Y0{x)^^^^Yl{x) + 
dx dx 

+ 4Y2{x)Y0{x)a2{x)Yl{x)Xl{x) - AY2{x)Y0{x)Yl{x)a2{x) '^^^'^^^ _ 

dx 

- 4Y2{x)Yl(xfa2{x)X0{x) - ma2{x)X4{x)Y0{xf + 

+ 8a2{x)X3{x)Yl{x)Y0{xf ~ 4Y^{x)a2{x)Y\{xf X2{x) - 

~ 2y0(x)ri(a;)3^^^^ + a2{x)Y\{xf X\{x)) j 
dx 

[\^Yf^i(xfX4{x)Y2(x) - 4Y<^{x)X4(x)Y\(xf), (43) 
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bO{x) = {Yl{xfa2{x)b2{x)^^^^ - 4Y2{x)Y0{x)Yl{x)a2{x)b2{x)^:^^^ + 

+ SY2{x)YG{xfb2(x) '^"''^y^ Yl(x) - 8Y2(x)Yl(x)Y0(xfX4(x) + 
dx 

+ AY2{x)YQ{x)a2{x)Yl{x)h2{x)Xl{x) - AYQ{x)h2{x)a2{x)Yl{xf X2{x) - 

- AY2{x)Yl{xfh2{x)a2{x)XQ{x) - \&a2{x)h2{x)X'i{x)YQ(xf + 
+ ?>a2{x)h2{x)Xi{x)Yl{x)YQ(xf + 2YQ{x)Yl{xfXA{x) - 

- 2rO(,T)fe2(x)ri(x)3^4^ + b2{x)a2{x)Yl{xfXl{x))l 

dx 

{-'iY{){x)a2{x)XA{x)Yl{xf + lQY{){xfa2{x)X'i{x)Y2{x)), (44) 



d'^a2{x) 



dx^ 

= {8a2{x)Y0{x)XA{x)X0{x)X3{x)Yl{x) +8Y2{x)a2{x)Y0{x)^Xl{x) 



dX4{x) 
dx 

- 2a2{x)Y0{x)XA{x)Yl{xf^^^^ - 8a2{x)Y0{xfYl{x)X2{x)^^^^ + 

+ 8a2{x)Y0{xfYl{x)X4{x)^^^ + 2a2ix)Y0{x)Yl{xf^^^^^^ + 

dx dx dx 

+ 2a2(x)Y0(x)Yl(xf^^^^Xl(x)+4Y2(x)a2(x)Y0{x)Xl{xfXA(x) - 
dx 

- 16a2{x)Y0{xfX4{xfX0{x) - 16a2{x)Y0{xfX4{x)^^^^ - 

- a2(x)X4(x)Yl(xfXl(xf + a2ix)X4(x)i^^^^^^'fYl(xf - 

dx 

- 16X'i{x)X0{xfa2{x)Y2{xf - 16Y2{x)a2{x)Y0{x)X4:{x)X0{x)X2{x) + 
+ 8Y2{x)a2{x)X4{x)Xl{x)X0{x)Yl{x) + 16Y2{x)YQ{xf 



^da2{x) dX4{x) 



dx dx 

- AYOixfYlixf^-^^-^ + 16a2ix)Y0{xfX3{x)^^ - 

dx dx dx 



8Y2{x)a2{x)Y0{xyX4(x) — + 8Y2{x)a2{x)Y0{x)Yl{x)X4:{x 



dx dx 

- 4Y2{x)a2{x)Y0{x)X4{x){^:^^f - 8Y2{x)a2{x)Y0{xf^^^^^^^^ 

dx dx dx 

- 2a2(a:)yO(x)X4(x)yi(x)2^^^^ + 8F2(x)a2(x)yO(x)'X4(x) '^^^^^''^ 



dx^ dx"^ 

-8Y2{x)a2{x)Y0(x)Yl{x)X0{x)^^^)/ 

dx 

{-4Yl{xfYQ{xfXA{x) + 16Y0{xfXA{x)Y2{x)), (46) 
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dY2(x) 

= (-Yl(xfY0(x)X3(x) - Yl(xfXl(x)Y2(x) - 

ax 

- Yl(xf^^^^Y2(x) + 4Yl(x)Y0(x)Y2(x)X2(x) + AYl(x)YO(x)^ XAix) + 
ax 

+ 'iYl{x)Y2{xfX0{x) - 'iY0{x)Y2{x)^Xl{x) + 

+ 4rO(a;)y2(x)2^^I^ - 4Y2(x)Y0(xfX3(x))/ 
ax 

{4:Y0{xfY2{x) - Yl{x)^YO{x)), (47) 

^IIM. = {8X0{x)Y0{x)Y2{x)^ - AY2{x)^^^^^Yl{x)Y0{x) - 
ax ax 

- 4Y2{x)X0{x)Yl{xf - 8X4:{x)Y0{xf + 4:X3{x)Yl{x)Y0{xf - 

- 2Y0{x)Yl{xfX2ix) + Yl{xfXl{x) + Ylixf^^^^^^)/ 
{~4Y0{xfY2{x) + Yl{x)^YO{x))l (48) 

where [a2{x) ^ 0, X4:{x) ^ 0, YO{x) ^ 0, 4Y0{x)Y2{x) - Yl{xf ^ 0] (49) 

The equations (^7|), (^ ) arc intcgrability conditions for this case. But equa- 
tions for undetermined parameters of the structure ( p7|) include ubiquitous linear 
second-order equation (|4^). From restriction ( ^9|) we see that a2{x) ^ 0, so the 
particular solution a2{x) = is not suitable in this case and in fact we have to 
solve (^6|) by existent methods. It is obvious that under some conditions the 
solutions of the equation (Efl) may involve special functions. 



7 Conclusions 

We presented some aspects in solving rational first-order ODE's by examination 
of the structures of general solutions. Examples in sections 4-6 demonstrated 
that proposed method gives as a first step the systematic way for obtaining new 
integrability conditions and general solutions for various families of rational 
first-order ODE's. 

Besides first-order ODE the method is suited for solving non-linear equations 
of higher order too. The properties of its solutions expect further studies. 
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